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o : 

^N) , For a non-interacting Bose gas on a lattice we compute the shift of the critical temperature 

for condensation when random-bond and onsite disorder are present. We evidence that the shift 
depends on the space dimensionality D and the filling fraction /. For D — > oo (infinite-range model), 
h-^ , using results from the theory of random matrices, we show that the shift of the critical temperature 

. is negative, depends on /, and vanishes only for large /. The connections with analogous results 

QQ ' obtained for the spherical model are discussed. For D = 3 we find that, for large /, the critical 

^ temperature Tc is enhanced by disorder and that the relative shift does not sensibly depend on /; 

at variance, for small /, Tc decreases in agreement with the results obtained for a Bose gas in the 
continuum. We also provide numerical estimates for the shift of the critical temperature due to 
' disorder induced on a non-interacting Bose gas by a bichromatic incommensurate potential. 
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I. INTRODUCTION 



The ideal Bose gas (IBG) is a well studied topic in quantum statistical mechanics [H, 0,111 since it provides invaluable 
qualitative insights on the properties of realistic bosonic systems. Its thermodynamics properties can be explicitly 
worked out and have been intensively studied both in the thermodynamical limit and in finite systems Nowadays, 
the study of the IBG has more than just academic interest since recent progresses in manipulating cold atomic 
I , systems 0, @ allow for a reliable control of the relative strength of the kinetic and interaction energies: for instance, 
■ the possibilities of varying the scattering length through Feshbach resonances [1, @ and superimposing optical lattices 
^ ' (obtaining effective lattice Hamiltonians 0,11]) provide a tool for testing a Bose gas from the strongly interacting 
O regime to the (almost) non-interacting one. In Ref. it has been shown that it ispossible to tune the s-wave 
scattering length a of a '^^Rb Bose gas with high precision, and set it almost to zero flOl : this allows for long-lived 
Bloch oscillations in an optical lattice [ll| and allows one to perform atom interferometry [ll| . When an optical lattice 
04 ■ is superimposed to an IBG, for large values of the laser power, the system is described by a tight-binding bosonic 
K*" ' model. 

I The IBG has also been used to understand the effect of disorder [ll] and the role played by confining, tailorable 
' external potentials, relevant in investigations of trapped ultracold atoms [1, The effects of controllable or random 
impurities is again motivated by the huge progresses made on the manipulation of ultracold bosonic gases culminated 
in the recent experimental observation of Anderson localization for matter waves in a random potential (l^ . ITsj . In a 
Bose gas, disorder can be induced both by a laser speckle (i-e., a laser beam is shone through a diffusive plate 

and superimposed to the trapping potential) or by an incommensurate bichromatic potential (i-e., an auxiliary 
, lattice is added to the main lattice, with an incommensurate ratio between the two frequencies). Random on-site 
disorder in the strongly interacting limit induces a Bose glass phase [H, [1^, whose characterization has been the 
object of many theoretical [2l|,[2l,[2l,[23,[2l,[2l,[23,[H,[H,l3l] and experimental [ll] investigations. Finally, ultracold 
atoms are a good experimental setup for studying sup erfluidity in random environments, which is an issue related to 
the adsorption of ^He in porous media [H, [U [H, [H, [H, HI, Isl] . 

The IBG is very interesting also in relation with statistical mechanics models: indeed, the IBG is in the universality 
class [37] of the spherical model (SM) [s^, which displays a phase transition in dimensions larger or equal than 
three. The SM is solvable in any dimension also with an external magnetic field js^Jsol \id\, and it provides an 



interesting laboratory to investigate a variety of topics, such as long-range interactions [41l | and finite size scaling 42 1. 
Furthermore, it is the n — > oo limit of the 0{n) model [ill, and, thus, is the leading term of 1 /n expansions 44 1. 
With disorder, the spin-glass SM with infinite-range interactions has been studied intensively [H, \^ and it has been 
used to analyze the properties of relaxor ferroelectrics 47] . 

When there is disorder, non-interacting bosons on a lattice are described by the Hamiltonian 



2 



with 



^ij — + Cij)- 



(2) 



In Eqs. ((I])-(l2l), the lattice sites are denoted by i,j, t is the tunnehng rate between neighboring sites, Aij is the 
adjacency matrix characterizing the lattice (i.e., Aij = 1 if i and j are nearest neighbour sites, and if they are not) 
and bi, h\ are the bosonic operators. The number of sites on a cubic lattice of dimension D is Ns ~ , while the 
total number of particles is Nt- The filling (i.e. the average number of particles per site) is given by 



Random on-site disorder is accounted by e^, while random-bond disorder is included in the definition of Aij through 
Cij; in the following Ci and Cy are regarded as random variables with vanishing average and variance v^t^ and v"^ 
respectively. We restrict ourself to situations where Uj = tAij > (u^ ^1), since they naturally arise when bosonic 
lattices are disordered by inserting a speckle or by adding an incommensurate optical lattice: in these situations both 
random-bond disorder Cij and random on-site energies occur and, typically, Vp > f. The inclusion of random-bond 
disorder allows also to make contact interesting properties of the spin-glass SM [4^ . 

In this paper we investigate how the critical temperature of a lattice IBG is shifted by the presence of disorder, 
elucidating the dependence of this shift on the filling. For the continuous Bose gas, the shift in the critical temperature 
has been derived in Ref. [sF]. In that paper it is shown that the critical temperature decreases with disorder: as we 
shall see, the lattice computation reproduces this behaviour for small filling, but, on the lattice, different behaviours 
are possible as a result of the interplay between discreteness and disorder. Indeed, on a lattice, there are two competing 
effects since disorder broadens the distribution of the delocalized eigenstates of Hamiltonian ([1]) and further localizes 
the ground-state wavefunction: the former (latter) effect tends to lower (increase) the critical temperature. For 
D ~* oo (i.e., infinite-range model), the analysis can be carried analytically and one finds that the two effects exactly 
balance at large filling while, for small filling, the critical temperature decreases. For D = 3 the balance is obtained 
for intermediate filling while, at large (small) filling, the critical temperature increases (decreases); the numerical 
coefficient that, for small disorder and large filling, determines the critical temperature increase is much smaller for 
random on-site disorder. 

The plan of the paper is the following: in Section II we determine the critical temperature of non-interacting 
bosons on a three-dimensional lattice; there, we analyze the limit of large filling and the connection with the SM. 
Section III is devoted to the analysis of the infinite-range model without and with disorder. In Section IV we study 
three-dimensional lattices with random-bond disorder and random on-site energies; in addition, we exhibit numerical 
results for the shift of the critical temperature due to disorder induced by a bichromatic incommensurate potential. 
Our concluding remarks are in Section V. 



In this Section we determine the critical temperature for condensation of bosons hopping without disorder on a 
three dimensional lattice. We shall show also that, at large filling, the value of this temperature coincides with the 
one of a mean spherical model. 



In the absence of disorder (e^ = and Cy — 0, Vi, j), the Hamiltonian ([T]) reduces to the boson- hopping Hamiltonian 



II. NON-INTERACTING BOSONS ON A LATTICE 



A. Hopping Bosons without Disorder 



Hio) 




(4) 



By setting 4 = Ng^^'^ J^j hjC^'^ , H^"\ H^^^ may be readily diagonalized; iov D ~ i one gets 




(5) 



k 



with 



2t (cos kx + cos ky -t- cos kz) ■ 



(6) 
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We assumed periodic boundary conditions, so that kg = 2TTne/L with = 0, • • • ,L — 1 {£ = x,y, z). 

In computing the grand-canonical partition function, one adds the chemical potential /i enforcing the conservation 
of the total number of particles; as a result, one replaces -ff with 7^'°) — H^'^^ — ^N, where N — J^i — dl^du 
is the operator counting the total number of particles. Requiring (N) = Nt yields 

E ,,iE.l) _ 1 - ^T, (7) 
k 

which fixes the chemical potential /i as a function of T and Nt- In the thermodynamical limit, Ns,Nt — > cxd at 
fixed fining /, Eq. ^ reads /g^ e^<^k->^)-i ~ '^^ere /g^ dk - ■ ■ denotes the integral over the first Brillouin 

zone: J^^ dk - ■ ■ — J^^ dk^dkydhz • • • . The critical temperature in the absence of disorder is denoted by Ti"-* and it 

is defined by the condition yu(Ti'^^) = — Eq, where Eq — —6t is the ground-state energy: 

/ (8) 



with (3i°^ = \/kBT^°\ 

A simple estimate of TjP^ may be obtained by keeping only the lowest order of the Taylor expansion of the exponential 
in Eq. One gets 



yielding 



In Eq. (Uni) W(l) ~ 1.516386 and 



^ ^—^f (9) 



fc„T(°) ~ (10) 

^ "w(i)- ^ ' 



W(z)= / \ (11) 

Jbz (2^)^ l-^EtiCOsA:, 



is the three-dimensional generalized Watson's integral [40]. Some useful properties of the generalized Watson's integral 
are grouped in Appendix A. 

In order to show that Eq. (jlOp provides indeed an accurate estimate of Ti*^^ for large filling, it is most convenient 
to write the exponential in Eq. ([5]) as a geometric series; this yields 

00 

j = ^e-6,Xo^^^(2jXo)f , (12) 
where Xq = and Iq is the modified Bessel function of the first kind [i^. When Xq <C 1, one can approximate 



the series in (IT^ with an integral and - for J 1 - one obtains / ~ (l/6A"o) /q dxe~^ [/o(x/3)] . From Eq. (|A2p . 
one readily gets Eq. (fTO|) . 

Eq. (jl2p is very useful to numerically compute Xq as a function of / (or vice versa). Indeed, although the series in 
(fT^ is slowly convergent, it is always possible to choose a large enough integer Af so that Eq. may be written as 

AT 

/ = E (2jXo)]' + n{Af), (13) 

where the remainder of the series TZ may be approximated (with arbitrary precision for a pertinent choice of Af) as 

In Eq. ([T4|) the Riemann zeta function C and the asymptotic expansion of the modified Bessel function of the first 
kind /o(a;) — /\/2ttx valid for large x [i^ have been used. 
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In Fig. [T] we plot the ratio Q between the critical temperature Tc^\ determined numerically from Eqs. ([13]) and 
(fTi)) . and the large filling critical temperature 6tf /kB\N{l) given by Eq. PH)) . as a function of the filling /. The solid 
line corresponds to the 1// expansion, which is found to be 



(0) 



6tf 
W(l) 



1 

27' 



(15) 



Notice that, at finite filling, Ti°^ is larger than the critical temperature given for large filling in Eq. ([T 
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FIG. 1: Open circles: ratio Q between the critical temperature ri"' and the large filling critical temperature 
function of the filling /. Solid line: expansion (|15p . 

We conclude this Section by observing that, using [ioj 



6t/ 
fcBW(l) 



lim — 



r^o- ^ — ' iw„ — X ef^^ — 1 

n— — oo 



with P = (3^c^ (the sum is over all the Matsubara frequencies w„ = 2-Knj 0^c^ and fi, = 1), Eq. ([HI) may be written as 

dk 1 



,(0) 



/ ^ kslf^ lim ^-^^"^ / # 

T^o- ^ Jbz (27r 



(16) 



which is valid for all values of /. From Eq. (|16p . Eq. ([9|) may be obtained by taking only the contribution from 
the lowest Matsubara frequency ujq. Contributions from higher Matsubara frequencies {n = ±1,±2, • • •) lead to a 

nonlinear equation relating ri"' and /. 

As we shall see, Eq. ([T0|) provides also the critical temperature of the SM discussed next. 



B. Correspondence with the Mean Spherical Model 



The Hamiltonian of the SM reads 



(17) 



where the values of the (classical) spins are real and satisfy the constraint J2i = ^s- In the original solution of the 
SM [sl] , the constraint is taken into account by writing the partition function as Z — j Y\.i dSie~^^ S{^^ Sf — Ns). At 
variance, in the so-called mean SM the constraint is satisfied only on average: upon introducing the grand-canonical 
Hamiltonian K = H — fJ,{J2i ^1 ~ ^s), one only requires {J^i Sf) = Ns, where (• • • ) denotes the thermal average 
with respect to K. 
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In order to compare the results of the mean SM with the ones of an IBG with definite filhng /, it is most convenient 
to write the constraint as: 

J^Sf^fNs^Nr. (18) 

i 

The grand-canonical Hamiltonian is 

K = H-^i(^Sf-fNsy (19) 

while the partition function is given by 

Z= f '[[dSje-'^^. (20) 
For a cubic three-dimensional lattice, performing the gaussian integration in Eq. (j20p yields 



k 

Minimizing the free energy F — —kBTlnZ with respect to /i leads to 



Z = e-f"'f^^'Y\ J—^ ^. (21) 



from which, in the thermodynamical limit, one gets 



The critical temperature Tc is given by the condition /i(T'c ) ~ Eo/2 yielding 

ksTi" = (24) 
Eq. (IM)) coincides with ^TU\i . which has been obtained for the lattice IBG for large filling. 

III. INFINITE-RANGE MODEL 

In this Section we investigate the D oo limit of the boson-hopping Hamiltonian ^ without and with random- 
bond disorder. The main advantage of the infinite-range model lies both on its solvability and on the insights it 
provides for the behaviour of a realistic three-dimensional lattice model. 

A. Without disorder 

Let us start by considering the infinite-range boson hopping model {Aij ~ I if i j) . Its Hamiltonian is given by 

The single-particle energies are straightforwardly computed: one finds a non-degenerate ground-state energy Eq — 
—t{l — 1/Ns) and an excited state Ei = t/Ns with degeneracy Ns — 1. Conservation of the total number of particles 
yields 

" e/3(Bo-p) _ I + e/3(-Bi-M) _ i' (2^^ 
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As usual, the critical temperature Tc'^^ is obtained from the condition fi{Tc^^) = Eq which, in the thermodynamical 
limit, yields 



fc^ri") = / (27) 



In ( 1 + ) 



For / > 1, Eq. ^ leads to 



fesTi") ~ tf. (28) 



A plot of the ratio Q between the critical temperature T^"'' [given by Eq. ((27|) ] and the large filling critical temper- 
ature [given by Eq. (|28|) ] is provided in Fig.O The corrections to the large filling limit are of order 1//; indeed 



fcBTi"^ _ 1 _ 1 

tf /ln(l + i) 



(29) 



As for the three-dimensional lattice, the critical temperature at finite filling given by Eq. ((27)) is larger than the critical 
temperature determined for large / in Eq. (|28|): it turns out that the 1// correction has the same coefficient as in 
Eq. (fT5|) . From Fig. [2] one sees that the expansion ((29|) is rather good also for intermediate filling, / > 1. 
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FIG. 2: Plot of the ratio Q between the critical temperature Ti and tf /ks as a function of the filling / (solid line). The 
dashed line corresponds to the expansion at order 1/f given by Eq. (|29p . 

To conclude this Section, we observe that if one adds a diagonal term to the Hamiltonian (P5)) {An — 1), one still 
finds that the critical temperature is determined by Eq. (|27p . 



B. With random-bond disorder 

In the presence of random-bond disorder the Hamiltonian is obtained from Eq. (|25p replacing Aij with Ai^ . Namely, 

^ = -j^ E(i + ^ - E J^Mh, (30) 

where, as in Ref. [45| . Jij is a real symmetric matrix with elements normally distributed, having average value t/Ns 
and variance J'^/Ns (but variance 2J'^/Ns on the diagonal). We set 

v^i (31) 

and we limit ourself to u < 1. The Hamiltonian (|30p provides a mean- field theory description of the bond-disordered 
boson hopping Hamiltonian discussed in Section IV. 
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In this Section we compute the dependence on the fiUing / of the shift of the critical temperature, 5Tc = Tc~ ri"-*, 
with Ti"^ defined in 1^. We find that while ST^ < for finite /, and, ST^ for / ^ oo. 

In the infinite-range model the density of states can be obtained through an orthogonal transformation mapping 
the problem in that of a single local impurity in a crystal (45l . |50| : one obtains the well-known semi-circular law j5l| 
and a single isolated impurity state which splits off from the continuous part of the spectrum. The same result can 
be obtained by the repHca trick The density of states p{E) is given by 



piE) = po{E) 



1 



5{E - So), 



where 



Eq 



The density of states po{E) is given by po{E) = for | i? |> 2J and 



(32) 



(33) 



(34) 



for \ E \< 2J. A plot of the spectrum is provided in Fig. [3l one sees that the ground-state energy decreases (this 
induces an increase of the critical temperature Tc for the bond-disordered model), while there is a spreading in the 
distribution of the continuous spectrum (instead, this induces a decrease of Tc). The shift in the critical temperature 
is a result of the competition of these two effects. 
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FIG. 3: Solid line: density of states p{E) (in units of l/vrj) as a function of E (in units of t) for the infinite-range model with 
V = 0.25 - the ground-state energy (|33p is represented by a vertical solid line. The eigenvalues E — and E — ~t in the absence 
of disorder (v — 0) are indicated by vertical dashed lines. 



Using Eq. (|32p one finds that the critical temperature is determined by (J3c ~ ^/kBTc): 

1 



/ = 



27rJ2 ./_2,7 e^=(^-^o) _ 1' 



(35) 



with Eq given by Eq. Defining X = (3ct and x = E/2J, Eq. becomes 

TT J .^ ^X{l+v^)+2Xvx _ I ^ 

In order to compute STc for / ^ f , one can, as in Section II, Taylor expand the exponential in the r.h.s. of Eq. ((36|) 
to get 



I 



dx- 



_ 1 

From Eq. (|28p . one immediately sees that 5Tc — 0. 



(37) 
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For any finite value of the filling / and for small v, STc may be computed from the series expansion of the integrand 
of Eq. ([M)) : one obtains 



1 



1 



(e^ - 1)2 



3^-1 2 ^ 



Upon defining = (3i°h = In (1 + 1//) and putting X = Xa + SX , one may readily linearize Eq. 
One finds 



SX ~ v^Xo 



Xoe 



Xo 
2 



from which one gets 



6Tr 



(0) 



(38) 
around Xq. 

(39) 
(40) 



Thus (JTc < 0. 

Eq. (1^0)) is plotted in Fig. U) In Fig. O we compare the result for 5Tc obtained from Eq. (UOl) with the numerical 
solution of Eq. ([55]) for different values of the filling and the disorder. Although Eq. is valid for small a rather 
good agreement is found also for relatively large values of the disorder parameter v. 

In order to compute / to any order in w, one observes that Eq. p6p can be written as 



p=i 



Ii{2vXp) ^_x(i+^2)j, 
2vXp 



(41) 



where Ii is a modified Bessel function of the first kind [i^. Expanding both Ii{2vXp) and e ^(^+'"^)p in terms of v 
and using the definition of polylogarithm function, Lin{x) — J2^=i f^; one gets 



fe,r!=0 



(-1)"X 



n!fc!r(fc + 2) 



-Li, 



(2fc+r 



,,2fc+2n 



(42) 



where F is the Gamma function. It is easy to verify that, up to the order v'^, Eq. ([42]) reduces to Eq. p8|) . 

We conclude this Section by noticing that, for the D oo SM, the effect of random-bond disorder has been worked 
out long ago [i^: for disorder smaller than a critical value, i.e. for w < 1, the authors of Ref. [iH found that the 
critical temperature Tc of the ferro-paramagnetic transition does not get shifted as a result of bond-disorder. Our 
result remarkably parallels the result obtained in Ref. [45| . 




FIG. 4: Shift of the critical temperature (in units of v'^Tc''^) as a function of the filling according Eq. (|40|) . tJ"' is the critical 
temperature for the lattice without bond-disorder and it is given by Eq. (|27p . 
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FIG. 5: Shift of the critical temperature 5Tc (in units of Tc ) as a function of the disorder v. circles (triangles) are obtained 
from the numerical solution of Eq. (I35|l for f — 1 {f — 0.1) - solid lines: analytical prediction Eq. (|40|) . Inset: STc/Tc'^^ from 
Eq. (|1D)) vs. the filling f ior v — 0.25 (squares) and v = 0.5 (filled circles) - solid lines: Eq. (gDJ. 



IV. THREE-DIMENSIONAL HOPPING BOSONS WITH BOND AND ON-SITE DISORDER 

In this Section we determine ST^ for non-interacting bosons m D — 3 with random-bond and on-site disorder. The 
Hamiltonian is given in Eq. ([1]), where Cy (ej) is a gaussian random variable with zero average and variance v'^ (w^t^). 
In the following we use the replica trick [46l [ssj : we shall determine the Green functions of the disordered system 
from the replicated action, where disorder enters as an effective attractive interaction between replicas. As usual, one 
introduces N replicas and, after averaging the partition function on the disorder, perform the limit iV ^ 0. In the 
following, for the sake of simplicity, we separately consider the effect of random-bond and on-site disorder. 



A. Random-Bond Disorder 



A standard analysis (reported in Appendix B) allows us, after integrating over disorder, to derive an effective replica 
action as 

^^ff = E r^ME^r(r)(|:-/i)^?(r)-iX^V'r(r)A,,^,"(r)i (43) 

a. ^ \ i ^ ij J 

-^E^^^-/^^ |''dr'(^rM¥',"(TVr(T0v'J(r') + ^r(T)^,"(r)^7(r')¥'7(T')). 

The fields ip^s represent Bose particles and a,j — 1, ...,N are the replica indices. As one can see from (|43|) . bond- 
disorder induces an effective interaction with the two-particle vertex functions depicted in Fig. [51 Eq. (I43p is derived 
under the assumption that Aij = Aji for any realization of disorder. 

a I y a / y 

-< 1 < ^ 1 <: 

I I 
I I 

> ^ ^ ^ 

a ./ y a ./ y 



FIG. 6: Two-particle vertices for the disorder induced effective interaction. 



10 



Defining the Green's functions 



Bz (27r)3 - z 



(44) 
(45) 



with z = /i + iujn and performing a perturbative evaluation of Eq. (|45D up to the first order in , one gets for N 



2j-2 



G, 



^zj Qjrn H~ Qjm Qii 



Upon introducing the self-energy function II(k, z), one may write the Dyson equation in Fourier space as follows 



where 



G-i(k,z) = g-i(k,z) + [Z(k,z)-,5M], 



g-i(k,z) = ^k-z 



and (5/^, the shift of the chemical potential, defined as 

= II(k = 0, z = ^). 

Bose-Einstein condensation occurs when 



G(k,z)-i 



k— 0,^— 



which is solved for 



-m- ^i+ [51(0, /i) - (^Ai] = 



(46) 

3 

(47) 
(48) 
(49) 
(50) 



thus, for bond-disorder, the chemical potential at Tc is given by /ic + 5^c ~ —dt + 51(0, ~6t). 

Upon Fourier transforming Eq. (j46[) , a straightforward but lengthy computation yields the Fourier component of Z 
up to the order w^. One gets 



Z(k,z) = -v^t^ 



1 



dk' X]i=i cos(fc,' - ki) 



(51) 



Elementary computations based on the definition of Watson integral (see Eq. (jlip ) show that, introducing the functions 

(52) 



S,{z)^ie I ^_^ = -^W(-z/6t), 
Jbz (27i-)3 



S2(^) 



t 



^3 i?k. - Z 

dk' Ei,, 



t 



2 (2^)^ E],,-z 2 
and using trigonometric relations, Eq. (j5ip may be compactly presented as 

1 



1- W(-z/6t) 



Z(k, z) = ( Si(z)- -S2(z)£;k 



(53) 



(54) 



From Eq. (j54p one may readily compute (5^c, as well as the renormalization of the bandwidth, which, at the critical 
point, are given by 



5^ic = -w^ (Si(-6<) + S2(-6t)) = w2(2W(l) - l)i/2 w w^t, 



tc = t + v^-S2{-6t) = t 
6 



l + _^;2(W(l)_l) 



(55) 
(56) 



Eq. ([55)1 shows that, as a result of the bond disorder, the bandwidth is increased: this is consistent with the mean-field 
result obtained in Section III. 
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The starting point of the computation of STc is given by Eq. (fT6|) . One notices that, as a result of the contribution 
in Eq. ([50)1 of the self-energy induced by the random-bond disorder, the filling fraction / is changed to / -I- Sf, with 
Sf given by 

dk ( 1 1 



Sf 



lim fcsTc 



— lim ksTc 



BZ 



dk [HCk, fic +iuJn) - SfXc] 
(27r)3 {Ei^ - - «w„)2 



Substituting 



dk 



BZ (2^)3 (£;k - z) 
dk 



(2^)3 (£;k - ^) 



1 A 

dz 



Si(z) 

S2W 



in Eq. ((57| one gets 

(5/ = lim fesTcy^ e' 



3^2 



Si(z) - Si(^c 



d 



—Si{z) + S2{z)—S2{z) 
oz oz 



d 
oz 



(57) 

(58) 
(59) 

(60) 



For large values of the filling /, the dominant contribution in Eq. (j60p is given by the lowest Matsubara frequency. 
The explicit computation carried in Appendix C yields 



Sf 



9 
4^ 



W(1)(W(1 



1 



-knTr 



-0.04— fcsTe, 



(61) 



As a result of Eq. (fTU]) . one notices that, for large /, there is a linear relation between ri°'' and /, which is given by 

W(l) 



6t 



■kBT^'^^f. 



At the order u^, taking into account the shift in the particle density given by Eq. ([61]l . one also has 

W(l) 



ksTc + Sf = f. 



From Eqs. ((62)) and ((63|) and using Eq. ((6T|) . one obtains 



Te = ri°)|i + J.2(^ 
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1 



implying that 



JTc/ri"' is then positive and independent on the filling, provided that / ^ 1. 

We numerically checked the validity of Eq. ([65]) . We adopted the following procedure. For several values of v and 
/ and for several (100) configurations of hopping disorder and finite sized samples (cubic lattice with L sites per side 
with L = 5, 10), we diagonalized exactly the Hamiltonian H{L), then we averaged the single-particle eigenvalues 
over all the configurations, computing Tc{L). Finally, we performed a finite-size scaling analysis using the equation 
YjTj ~ ~ T in order to fit the parameters c and = Tc(oo). 

In Fig. [7| the numerical results obtained for different values of v and / are compared with the analytic result (|65|) . 
One sees that, for small disorder and / > 1, there is a good agreement between the two results. For / very small, 
instead, we recover the result obtained by Lopatin and Vinokur 35]: indeed, for small filling, the lattice details become 
irrelevant and one may treat the problem as in the continuum limit. Using the results of 35] with t ^ /2raa'^ (ag 

AT' 2 

being the lattice constant), one has w — 0.05 The resulting curve for / = 0.01 is in very good agreement 

with the numerical results. The reason for this agreement is that, for a small filling, the occupied single-particle states 
lie at the bottom of the first energy band rendering the lattice ideal Bose gas basically equivalent to the continuous 
model. 

In Appendix D we consider also the effect of an asymmetric random-bond disorder (c.y 7^ Cji): the final result (|D4p 
is similar to (|65p . with a different numerical coefficient. 
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FIG. 7: Shift of the critical temperature STc (in units of the critical temperature of the ordered lattice) as a function of the 
disorder strength v. The dots are obtained by numerical analysis while the two dash- lines are the analytic results: respectively, 
from the top of the figure, our result (I65p valid for large / and the result of [S^ for the continuous disordered Bose gas. 



B. Effects of on-site disorder 



It is easy to generalize our results to situations where also on-site disorder is present, see Eq. ((T|). Introducing the 
parameter 7, determined by the delta-correlated on-site random energies Ci with variance v'^t'^, as 

_ Sv^ + vl 

7=^^, (66) 

one gets that Eq. is modified as 

Z(k, z) = (^7Si(z) - ^S2{z)Ei}j . (67) 

When only on-site disorder is present {v"^ 0), one has that 5I(k, z) = —{vl/3)Si{z): this signifies that there is no 
renormalization of the bandwidth induced by a small random on-site disorder. Eq. (|67p amount to state that the first 
term in Eq. (j60p should be multiplied by 7. As a result, when both random hopping and on-site disorder are present, 
the critical temperature is given by 

Tc - |l + + 1 (W(l) - 1)) } « Ti«) (1 + 0.16 + 0.01 vl) . (68) 

For comparable disorder strengths, Tc is more sensitive to random-bond disorder, since the numerical coefiicient of 
in Eq. ((68|) is smaller than the one associated to tj^. As a result, the dependence of Tc on a small on-site disorder 
is extremely weak. We checked numerically that, for v — 0, Vp S, 1 and large /, 5Tc is very small. At variance, when 
/ ^ 1, 5Tc is negative in agreement with the results of Ref. [sF 



C. Incommensurate potentials 



An experimental strategy employed to simulate the effect of disorder is based to the use a bichromatic incom- 
mensurate lattice 14, 18]. The ultracold atomic gas is subjected to a periodic potential of the form V [v) — 
Vmain (i") + Vgec (i")- While the main lattice potential creates the ordered lattice [e.g., an isotropic three-dimensional 
lattice corresponding to Vmain (r) = Vq (sin^ {kox) + sin^ {koy) + sin^ (koz))], the addition of a secondary lattice, with 
a spacing incommensurate with the one of the main lattice, leads to the emergence of a quasi-random potential 
which, in a finite sample, is equivalent, for all practical purposes, to a random potential. If the secondary poten- 
tial is one-dimensional, T^ec (r) = Vi sin^ {kix), then, for large values of the height Vq of the main lattice, one gets 
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the Bose-Hubbard Hamiltonian ([T]) with on-site energies of the form e,; = —vcos{2Trqix), where v depends on Vi 
{v ~ Vi/2), q — fci/fco s-iid i = {ixi'i'yi'iz), with the i^'s integers. In Ref. [l8|, the experimental values of Aq = 27r/fco 
and Ai = 27r/fci are SSOnm and 1076nm, respectively. Of course, if the secondary lattice is three-dimensional, one 
would have ei oc v (cos {2'Kqix) -\- cos {2-Kqiy) -|- cos (2'Kqiz)). 

In Fig. [8] we plot, for two values of q, the shift in Tc for a finite size system (L ~ 10) as a function of the disorder 
intensity w of a one-dimensional secondary lattice for two different values of the filling: we see that, for small v, 
STc > while, for intermediate v, STc may become negative. For small v, STc oc w^. 



E^ 



-0.4' ' ' ' L 

1 2 



V 



FIG. 8: Shift of the critical temperature (in units of the critical temperature Tc'^'' (L) of the ordered lattice having a size L) 
as a function of v for a bichromatic incommensurate lattice. Open circles (open squares): ,f ~ i {f = 10) for q = 0.206. 
Diamonds (stars): f = I {f ~ 10) for q = 830/1076. In both cases L = 10. Dotted lines (for q — 0.206) and dashed lines (for 
q — 830/1076) are just a guide for the eye. 



V. CONCLUSIONS 

We investigated how the critical temperature for condensation of an ideal Bose gas on a cubic lattice is shifted as a 
consequence of random-bond and on-site disorder. Our analysis evidences how the shift in the critical temperature, 
6Tc, depends on the filling fraction / and on the space dimension D. 

For the infinite-range model {D — > oo), using results from the theory of random matrices, we found the explicit 
dependence on / of ST^. STc is negative and vanishes for large /, in agreement with known results for the SM 45]. 

For D = 3, we wrote down the corrections to the Green functions due to disorder, finding that, for large /, STc/Tc^^ 
is positive and does not sensibly depend on /. At variance, for small /, 6Tc is negative and our results are in a 
very good agreement with the ones obtained in Ref. [35j for the continuous disordered Bose gas. The reason for this 
agreement is that, for small filling (i.e, much smaller than one particle per site), the occupied single-particle states lie 
at the bottom of the first energy band: this renders the lattice ideal Bose gas basically equivalent to the continuous 
model investigated in Ref. [35j • 

We presented numerical results for the shift of the critical temperature induced by the disorder generated by 
a bichromatic incommensurate potential in a realistic finite size system: we found that the critical temperature 
increases for small values of the height v of the secondary lattice (i.e., small disorder), and decreases for larger values 
of V. The dependence of the shift of the critical temperature upon v is STc . 

We observe that, in the replicated action, the disorder enters as a fictitious effective attractive interaction between 
replicas. From this point of view, some similarity emerges with the problem of determining the shift of the critical 
temperature in a Bose gas due to the presence of a repulsive interaction [H, [13, IH, [H^ (for a clear presentation 
of this issue and more references, see [60|): rather counterintuitively, the critical temperature increases for small 
repulsive interactions, while decreases for larger interaction. This should be compared with the results obtained 
for an homogeneous continuous Bose gas with weak disorder (ssl . Isol . [6l| . where Tc decreases with disorder, as a 
consequence of the inter-replicas attraction. For the ideal Bose gas we showed in this paper that, on a lattice, this 
result is modified: indeed, for small filling, only the bottom part of the energy spectrum is occupied and the effective 
attraction among replicas amounts to a decrease of the critical temperature, in agreement with [35| . At variance, for 
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large filling, all the single-particle energy band is occupied, including states with effective negative mass, leading to a 
net increase of T^. 

We have considered non-interacting bosons on lattices with disorder: adding a repulsive interaction among the 
bosons, and treating it within the same approach we used, one expects [60l| that Tc is unchanged. Going beyond the 
mean-field approach of this paper should unveil the interplay between the discreteness induced by the lattice and the 
interaction in a disordered environment. For this purpose, non-perturbative resummation techniques to compute the 
self-energy, such as the one discussed in Ref. 62], could be relevant to determine 5Tc- In the same way, it would be 
very interesting to study the effects of an harmonic confinement on the disordered lattice i6vl . 

To conclude, we mention that if one introduces suitable time-dependent optical lattices jUll^, one should be able 
to see the effects also of negative effective 's: for non- interacting bosons, this could provide an experimental testing 
ground for the three-dimensional spin glass spherical model. 
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APPENDIX A: SOME USEFUL PROPERTIES OF THE WATSON INTEGRALS 

The generalized Watson's integrals are defined by 



(27r) J-.^^^ l-^Ej=iCOsfc, 



in our paper we are concerned only with W(2;) = W(3, z). 
The Watson integral pT|) admits an integral representation 

/•oo 

W(z)=/ dxe-^'lhix/iz)]^ (A2) 
Jo 

in terms of the modified Bessel function of the first kind /q- Particularly important is its value at z = 1, which is 
given by W(l) = (i) F (A) r (^) F (ii) 1.516386 0. 

The Watson function W(z) can also be written as a probability generating function 

oo 

WW = E^' 1^1 ^1 (A3) 

n=0 

with Pn being the probability that a random walker returns after n steps at the same point of the cubic lattice. In 
particular, 

(3 \ 2" 

is the probability for a random walker on a cubic lattice of reaching the starting site after 2n steps. Clearly, P2n+i = 
since a closed loop on a square lattice has always an even number of segments. For this reason one can write W(— z/6i), 
which is needed in the computation of Section IV and Appendix C, as 

n=0 

We notice that W(z) is also the solution of a Riemann p-differential equation, so that, when \z\ > 6t, \N{—z/6t) and 
its derivative are well defined. 
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APPENDIX B: EFFECTIVE REPLICATED ACTION 

Let us consider the Hamiltonian 

H = -tJ2 ^^Mbj - M E ^1^^ + E 'Mk (Bl) 

i,j i i 

where Aij = Aij{l + Cij), with Aij = 1,0 depending on if i and j are neighbour sites or not, Cy, being random 
variables. 

The partition function can be written as a coherent-state path integral (in units where h = ks — 1) 

Z = J \{v^{D^*e-'^, (B2) 

i 

where S is given by 



h3 



The partition function (|B2[) depends on the disorder realization of the hopping {c} and of the on-site energies {e}: 
Z = Z{{c}, {e}). Labeling the N replicas by a = 1, • • • ,N and using Aij = Aji, one has 

If, for simplicity, one assumes that = (only random-bond disorder is present), the averaged effective partition 
function is 

W{{c}) = ll[dc,,P{c,,)Z''{{c}), (B5) 
where P{cij) is the gaussian probability distribution 



1 



2 /^,,2 



P(q,) = ^— , (B6) 

with variance i;^. After integrating over the c^ 's, one gets 

Wi{c})^ JUVipfVipre-^'", (B7) 

where Seff is given in Eq. (j43p . When on-site disorder {e} is present, one can follow the same procedure: namely, 
one averages Eq. (jB4[) over disorder with the measure JJideiPo{ei) = Yl^dei —^==6^^'^*^^/^''°*^, where v^i^ is the 

variance of on-site random energies {vo is a dimensionless quantity). 



APPENDIX C: SHIFT OF THE FILLING FRACTION 

Using Eqs. ^E\) and in order to write Sf only in terms of the Watson function W(-z/6t), one 

obtains Eq. ((60)) . In Eq. ((60|) the dominant contribution to Sf is given by the lowest Matsubara frequency (tj„ — 0): 

= ksT,^ hm I [Si(z) - Si(/.c)l |-Si(z) + S2(z) |-S2(z) - S2(Aic) |-Si(z) j . (CI) 
' J az oz oz \ 



3^2 



2: — >^c 



One can compute Eq. (jCip taking advantage of the asymptotic behavior of \N(z). For = z"^ — (6^)^ ^ one can 
write 



27r 



W(-z/6t) = W(l) - ^,/l-M! + ^ fl - ^ ) + 0(,3) (C2) 
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where /C = ^ (w(l) + ^^5W(T)) • From Eqs. one has 



Si(z) = — I W(l)-^Vl-^ )' (C3) 



S2{z) = -lh- W(l) + ^^1-^ ) + Oiv'). (C4) 



The derivative of these functions are 



^SA.)^,eh^(J^^^^}m_&£^^Oi,,), (C5) 

dz I 27r z4 L (6t)^ ^2 z"! ' 



One should notice that Eqs. (|C5[ IC6|) diverge at z = —6t, but these divergences are canceled in Eq. (jCip . Inserting 
Eqs. (IC4|-(IC6| in Eq. jCl]), one readily finds Eq. (jCT]) . 

An alternative way to obtain Eq. ((6T|) is based on the random-walk probability expansion of the Watson integral 
presented in Appendix A. If one truncates the sum (IA5|) . the result for STc would be zero since, only for the infinite 
sum, one has that 

In-.m^-t^t^ (C7) 

71—1 

is not convergent at z = —6t. Eq. (jC7|) implies that only very long random walks are important: thus, expansion in 
terms of n^^ (the inverse length of the walks) is meaningless since only walks with n —^ oo are relevant, although P„ 
decreases very fast with n. For large n, P2n has the following asymptotic behavior 

/ 3 

P.„^2(— j ,.^oo. (C8) 

Inserting Eq. ([CS]) in Eq. ((J5|l (for n > 1, and Po = 1) and using Eqs (|52[ (j53p and (|C1[) ). after taking the limit 
z — > — 6i, one gets Eq. (f6T|) . 

APPENDIX D: ASYMMETRIC RANDOM-BOND DISORDER 

Here we consider the situation where Cij and Cji are two independent random variables, namely 

c,j ^ Cji. (Dl) 

This situation may arise, for instance, when time reversal symmetry is broken. 
After integrating over disorder, one gets 

GlTu = Glm + vlt" J2 "^V (N + 1)^H Grj Gjm + 0{v\) (D2) 

(where is the variance of the asymmetric random-bond disorder), which implies that 

Z(k,z) = ^S2(z)i?k. (D3) 
As a result, an asymmetric random-bond disorder yields a value of Tc given by 

Te = {l + + 1 (w(l) - l)) } . r(") (1 + 0.244) . (D4) 



17 



In order to convince oneself that time reversal symmetry could play a role, one may consider a system with only 
two sites and random hopping between them, described by a matrix 



ii2 

t21 



(D5) 



The eigenvalues of (jD5|) are 



A± = ±Vh2t2^. (D6) 

If time reversal symmetry holds, one has that 

tl2 = t21. (D7) 
Averaging on the disorder and assuming ti2 = t2i = t, one gets 

Al=±t (D8) 
On the other hand, if time reversal symmetry is broken (i.e., ti2 ^ ^21)1 one has 



3e J \ 12, 

In deriving Eq. (|D9p a uniform distribution for ti2 and t2i with width 2te and mean value t has been assumed 
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